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Abstract 

The occurrence of swarms of small bubbles in a variety of industrial systems enhances their 
performance. However, the effects that size polydispersity may produce on the stability of kinematic 
waves, the gain factor, mean bubble velocity, kinematic and dynamic wave velocities is, to our 
knowledge, not yet well established. We found that size polydispersity enhances the stability of a 
bubble column by a factor of about 23% as a function of frequency and for a particular type of 
bubble column. In this way our model predicts effects that might be verified experimentally but 
this, however, remain to be assessed. Our results reinforce the point of view advocated in this 
work in the sense that a description of a bubble column based on the concept of randomness of a 
bubble cloud and average properties of the fluid motion, may be a useful approach that has not 
been exploited in engineering systems. 
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I. INTRODUCTION 



The theoretical description of multiphase flows is essentially based on analyzing the re- 
ponse of a cloud of dispersed particles of different size ranges in a fluid. These particles 
constitute dynamic phases and hence a multiphase flow. A widely used multiphase system 
is a bubble column which is a reactor where a discontinuous gas phase in the form of bubbles, 
moves relative to a continuous phase. Bubble columns have a wide range of applications in 
conical industries, biotechnology 0r in nuclear reactors Q, Q, Q, J g The transient 
behavior is important at the start-up of these systems and its analysis is essential in order 
to characterize the dynamic performance of the columns. Among the phenomena that occur 
in these systems void wave propagation mechanisms are of great importance since many 
transient and steady states are controlled by the propagation of these waves and, in this 
sense among others, the dynamic characterization of multiphase flows is essential for the 
prevention of instabilities. 

The (in) stability of bubbly flows which are characterized by almost uniformly sized bub- 
bles, is usually described in terms of the propagation properties of void fraction and pressure 
disturbances caused by natural or imposed fluctuations of the rate of air supply 
Bubble size, rise velocity, size distribution and liquid and bubble velocity profile have a direct 
bearing on the performance of bubble columns. However, most of the time the dispersion 
devices deliver a dispersed phase with a given size distribution. The importance of the size 
distributions is only scarcely evaluated, most of the time by direct empirical trials and its 
influence on the global behavior has still to be studied. Actually, to our knowledge and from 
the theoretical point of view, it has not been yet well established whether the stability of the 
motion of a swarm of bubbles is different for monodispersed or polydispersed bubble flows. 
The main objective of this work is to investigate the effects that size polydispersity might 
produce on the stability of a bubble column. We shall introduce the effect of polidispersity 
through the drag force in the hydrodynamic equations, using a method based on statistical 
concepts and on a point-force approximation As we shall see below, the corrections on 
the drag force factor, C£(a), due to polydispersity depend only on the first three moments 
of a given particle size distribution and they also have an effect on several properties of 
kinematic waves. In particular, we found that size polydispersity enhances the stability of 
void waves by a factor which varies between 4.5 — 23% as a function of frequency and for 
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a particular type of bubble column. In this way our model predicts effects that might be 
verified experimentally but this, however, remain to be assessed. 

To this end the paper is organized as follows. In Section 2 we briefly review a hydro- 
dynamical model for bubbly fluids introduced by Biesheuvel and Gorissen 10j. Next, in 
Section 3 we consider a dispersion of spherical air bubbles of different radi in water and we 
calculate the effect of size polidispersity on the gain factor, mean bubble velocity, kinematic 
wave velocities as a function of void fraction, for different wave frequencies. 



II. EQUATIONS OF MOTION OF A BUBBLE DISPERSION 



In this section we summarize the main ideas and steps behind the hydrodynamical model 
for bubbly fluids introduced in Ref. Q]. The equations of motion for a swarm of bubbles in 
a bubble column have been derived in the literature by using standard methods of kinetic 
theory to average over an ensamble or realizations of the flow Q, Q. In Ref. Q a 
dispersion of equally sized air bubbles in a water column where the bubbles are small enough 
to remain spherical through the whole system, is considered. They assumed that the air can 
be taken as an incompressible fluid where no mass transfer is allowed between the bubbles 
and the water, which is assumed to be an incompressible Newtonian liquid. The conservation 
equation for the mean number density of the gas bubbles, n, and the conservation equation 
for the mean bubble momentum, p G ~v (Kelvin impulse), were obtained for this system [r^ . 

ft. _ . . . 



dt 



+ V, 



mv 







d_ 
di 



» [ ^a 3 p G v + l L 



n ( -ira 3 p G v + 1 L 



nF D + n-Tra 3 (p L - p G )\ 



(2) 



Ix is the fluid impulse, Tl (x, t) and (x, t) are the fluid stresses; is the drag force 
exerted by the fluid on the bubble and g stands for the gravity field; pi, p G denote, respec- 
tively, the mass densities of water and air. pi stands for the liquid's viscosity. In order to 
describe the flow parameters of the bubble swarm, Eqs. (JIJ and (J2J) should be expressed 
in terms of the volume fraction of bubbles (or void fraction) e and their velocity field v. 
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Following Ref. 10] we assume that the uniform flow of bubbles is along the axial direction 
of the column with a mean axial rise velocity vq(e), Therefore, e(z,t) = ^7ra 3 n(z,t). 

The effect of hydro dynamic interactions between the bubbles on the mean frictional force 
may be represented by introducing a function fo(e) into vq(e) in the form vq(e) = (e^oc, 
The magnitude of the terminal velocity, t>oo, of a single bubble of radius a in a stagnant 
liquid is given by [r| Voo = C^ 1 {pl — Pg)<?> where Co = is the drag force factor and 

experiments suggest that [14], fo(e) = (1 — e)~ 2 ■ The mean fluid impulse is modelled by 
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nI L = n 0vra 3 pL^ m (e)v (e), (3) 
where mn (e) takes into account the effect of the hydrodynamic interactions. According to 



Ref. |15( an expression for m (e) that renders reliable results up to large values of e is 
mo(e) = (l + 2e)/(l-e). 

Since in a nonuniform bubbly flow the stress T = Tq+ 1l play the role of an effective 
pressure, they also assume that the kinetic contribution, p e (e), is proportional to the effective 
density of the bubbles, e _1 p e /(e) = pc + |pl^o(^), and to the mean square of their velocity 
fluctuations Av 2 = H{e)vl{e) = j- (l — -f-^j Vq(e), [lq. Here e cp stands for the limit of 
closest packaging of a set of spheres and is close to the value 0.62. Thus, p e (s) = pefAv 2 . 
Furthermore, if the non-uniformity is the main cause of an additional transfer of bubble 
momentum and fluid impulse associated with stress, Biesheuvel and Gorissen [3] postulate 
that such a contribution to the stress should be given by the force p e (e)^v. Therefore, 
taking into account both contributions to the stress, T = — p e (s) + /i e (e)|^, where v is the 
one dimensional nonuniform flow velocity and p e (e) = ap e f(e)v (e)H 1 ^ 2 (e) is an effective 
viscosity. 

On the other hand, the mean frictional force is enhanced by an effective diffusive flux 
of bubbles due to their fluctuating motion. This effect is similar to an steady drag force 
acting upon each one of the bubbles and proportional to the mean number density gradient. 
Therefore, using (J2J) this force is represented by nF D = C D efo(s)[v + ^§f]. Substitution 
of the above expressions into Eqs. (0) and (0) leads to the following closed set of one- 
dimensional equations of motion for the bubbly flow in a zero volume flux reference frame, 

« + &(»)=». w 
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9 [Pef(e)v] + ^ [ P ef{e)v 2 ] - J^T 



dt 



ep ef dz 



(5) 



These equations may be rewritten in a laboratory reference frame by considering the 
mean axial velocity of the dispersion, U, defined by U(t) = eUg + (1 — s) Ul- Here Uq and 
Ul are the mean bubble and fluid axial velocity in the laboratory reference frame. Note that 
due to the incompressibility of both, liquid and gas, U is only a function of time. Therefore 
v = Ug — U and a Galileo transformation of Eqs. (p£J) and (JHJ) gives 



de d 

dt + ir/ UG = > 



(6) 



d_ 

di 



£ [ PgU g + ^pLm (U G - U) 



d_ 

dz 



£ \ PgU g + ^p L m (U G - U) ) U g 



d ( ^ dU G 
dz \ dz 



~C D ef (U G -U) + 



dU 
~ £PG ~dt 
p e (e) de 
epef dz 



£ (pa ~ Pl) g, 



together with the incompressibility condition 



— 

dz 



(7) 



(8) 



Consider a quiscent equilibrium state of the dispersion described by e = e . The devia- 
tions from this state will be denoted by 5e(z,t) and Sv(z,t). Linearization of Eqs. (JHJ) - (jSJ) 
around the reference state yields the wave-hierarchy equation 



T, 
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dz 
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dt 
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1 ' dz 
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dz 
d 2 8e 
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d_ 
dt 



d \ d 2 5s 



J G - 



dz J dz 2 



(9) 



with lower and higher-order wave velocities given by c = Uq + Sqv' an d 



+ 



U 



Go 



-^eoPLVorriQ 
Pg + \pLm 



± 



\e p L v m^\ 2 



.Pg + Iplitiq J Pg + ^Pl^o 

Here v e (e) = v (e)H 1 ^ 2 (e) and r e (e) = (Cofo)~ 1 [pg + |pl^o(^)]- The primes (/) denote 
derivatives with respect to e and evaluated at the unperturbed state e = e . 
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Pa 



1/2 



(10) 



6 



For relatively low radial frequencies the wave propagation is described by a linearized 
Burgers/Korteweg-de Vries equation 



d 







dt ~^ °° dz 



d 2 6c d^Ss 
8e « [r e (c + - c )(c - c~) + z/ £ ] — - + t £ v £ (U Gq - c ) 



dz 2 b tv "° uy <9,2 3 
with a solution e oc exp (7,2 — iut) where u is the frequency of the void wave and 
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10J 

7 ftj 

c 



z/ £ r e u; 2 (C/ G - c ) 



v e u? (U G - c ) 



c )(c - c ) + 5 e ] . 



;i2) 



In terms of these quantities the so called gain factor, Gf = exp [Re(7)u; 2 A,2], where Re(7) 
denotes the real part and Az the distance between two impedance probes in the experiments 
to measure Gf [18] . 



III. POLYDISPERSED DISPERSION 

The method developed by Tarn [9'] uses the concept of randomness of the bubble cloud 
and derives equations describing the average properties of the fluid motion. These averages 
are taken over a statistical ensemble of particle configurations. A slow viscous flow past a 
large collection of spheres of a given size distribution is considered to derive a particle drag 
formula free from empirical assumptions. The result essentially replaces the disturbance 
produced by a sphere in low Reynolds number flow, by that of a point force located at the 
centre of the sphere. The correction drag force factor is given by 

1 



C D — \Cd 
where 



1 + aa H — (era) 2 
3 



C D , (13) 



6nM 2 + r(6vrM 2 ) 2 + 12ttMi(1 - 3c)] 1/2 . s 

a = . (14) 

(l-3c) 1 J 

M n = J n(a)a n da are the moments of the size distribution n(a) and c = |7rM 3 . 

Since the terminal velocity of a bubble depends on Co, it is reasonable to assume that 
in the polydispersed case vo(e) should be replaced by Vq = A _1 wo. Substitution of this 
asumption into Eqs. (JBJ) - ©, carrying out the linearization procedure described in the last 
section and using the explicit expressions of (3 = {Ug , c , c ± , t £ }, one can show that these 
quantities scale as (3 P = {A _1 f/ Go , A~ 1 c , A~ 1 c ± , A _1 r e }. If these polydispersed quantities 
are substituted into Eq. (|T^|) . one obtains an expression for the polydispersed gain factor 
G p f = exp [Re{Y)u 2 Az}. 
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FIG. 1: Gain factors Gf, G p vs. e for waves with frequencies of 2, 2.5 and 3Hz. The liquid is 
stagnant and the parameter values are those given in section 4. 

IV. RESULTS 

To compare the monodisperesed and polydispersed results on the gain factor, mean bubble 
velocity, kinematic wave velocities as a function of void fraction, we used the the following 
parameter values for an air-water bubble column, Az = 20 cm, Vr = 1280 cm 3 , pc = 1.2046 
x 10 -3 gr/cm 3 , p L = 0.998 gr/cm 3 , \i L = 1.002 x 10~ 2 poise, e cp = 0.62. In Fig. 1 we 
plot Gf and G p vs. e for different frequencies u and for a log- normal distribution n(a) with 
average and dispersion a = 0.04 cm, o = 0.5, respectively. 

Note that for values 0.185 ^ e ^ 0.301 the attenuation rate drops significantly for For 
instance,the per cent difference defined by Y = \Gf — G P j\ /Gf. For the range both cases. For 
instance, for a frequency of 2Hz this difference ranges from 0.1 — 4.98 per cent, whereas for 
a frequency of 3Hz it varies in the interval 0.1 — 22.78 per cent. This means that stability is 
larg er in about 23 per cent for the latter case, a change that is significant in bubble reactors 



The quantities (3 = {Ug , Co, c 1 * 1 , r e } and f3 p = {A -1 ^,,, A _1 co, A~ 1 c ± } are plotted as 
functions of e in Fig. 2. The curve for Co is always between that for c +p and c~ p . According 
to the Whitham stability criterion Q], when Co < c™ the uniform flow is unstable. This 
occurs for both distributions, however, for the monodispersed case it occurs for e > 0.353, 
whereas for the polydispersed case the system is stable up to a larger value of the void 
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FIG. 2: (3 and (5 P as functions of the void fraction e for the same parameter values as in Fig. 1. 
fraction, e.g. e > 0.36. 



V. DISCUSSION 

Summarizing, in this work we have analyzed the effects of size polydispersity in several 
features of the void fraction waves and their stability properties. We found that the presence 
of a size distribution reinforces the stability of the waves, as shown in Figs. 1 and 2. 
Furthermore, the per cent difference may be quantified by estimating T = \(3 — (3 p \//3, T 
amounts to a maximum percentual difference of 4.9%. 

It is convenient to emphasize once again, that the hydrodynamic model used in this 
work is idealized in many aspects. For instance, compressibility and hydrodynamic 
interactions between bubbles and with the boundaries, have not been taken into account. 
However, given the complexity of these effects and of the system itself, the simple dimensional 
model proposed by Biesheuvel and Gorissen seems to be a good first step in modeling the 
complex behavior of a bubble column. It also ilustrates how some of the methodology 
and concepts of kinetic theory and statistical mechanics may be used to deal with complex 
phenomena in engineering systems. 

We should also mention that in this work we have assumed an initial polydisperse size 
distribution and the coalescence of bubbles has not been considered [20] . however, this 
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remains to be assessed. Some other important effects remain to be considered as well, 
like the bubble-bubble interaction mechanisms. Nevertheless the approach followed here by 
including the influence of the distribution through the drag effects, considering a mean field 
approaches an attempt to set a first framework for the bubble size distribution incorporation 
to further studies. Our results reinforce this point of view in the sense that a description of a 
bubble column based on the concept of randomness of a bubble cloud and average properties 
of the fluid motion, may be a useful approach that has not been exploited in engineering 
systems. 
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